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Abstract 

The TT+TT^ s-wave scattering phase-shift is determined below the inelastic threshold using Lattice 
QCD. Calculations were performed at a pion mass of rnvr ~ 390 MeV with an anisotropic nj = 2-1-1 
clover fermion discretization in four lattice volumes, with spatial extent L ~ 2.0,2.5,3.0 and 
3.9 fm, and with a lattice spacing of 5s ~ 0.123 fm in the spatial direction and bt ~ bs/3.5 in 
the time direction. The phase-shift is determined from the energy-eigenvalues of tt~^tt^ systems 
with both zero and non-zero total momentum in the lattice volume using Liischer's method. Our 
calculations are precise enough to allow for a determination of the threshold scattering parameters, 
the scattering length a, the effective range r, and the shape-parameter P, in this channel and to 
examine the prediction of two-flavor chiral perturbation theory: rn^ar = 3 -|- ©(m^/A?,). Chiral 
perturbation theory is used, with the Lattice QCD results as input, to predict the scattering 
phase-shift (and threshold parameters) at the physical pion mass. Our results are consistent with 
determinations from the Roy equations and with the existing experimental phase shift data. 
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I. INTRODUCTION 

Pion-pion (tttt) scattering at low energies is the theoretically simplest and best-understood 
hadronic scattering process. Its simplicity and tractability follow from the pseudo-Goldstone 
boson nature of the pion, a consequence of the spontaneously broken chiral symmetry of 
QCD, which implies powerful constraints on its low-momentum interactions. The ampli- 
tudes for 7171 scattering are uniquely predicted at leading order (LO) in chiral perturbation 
theory (xPT) pLj. Subleading orders in the chiral expansion give rise to perturbatively-small 
deviations from the LO determinations (for small pion masses), and contain both calculable 
non-analytic contributions and analytic terms with low-energy constants (LEC's) that can- 
not be determined by chiral symmetry alone [2H1]. Fortunately, Lattice QCD calculations 
are reaching a level of precision where statistically significant values of the LEC's in the 
1 = 2 (7r"'"7r"'") channel are being calculated. Once the LEC's are obtained using unphysical 
lattice pion masses, xPT can be used to predict the phase shift at the physical pion mass to 
relatively high precision and with quantified uncertainties. The current capability of Lattice 
QCD — in conjuction with xPT — to calculate tttt scattering parameters very accurately is 
important theoretically because Roy-equation [SHTj determinations of tttt scattering param- 
eters, which use dispersion theory to relate scattering data at high energies to the scattering 
amplitude near threshold, have also reached a remarkable level of precision [814T0]. and the 
results of the two methods can now be compared and contrasted. 

There have been independent lattice QCD determinations of the 7r"^7r+ scattering length; 
with three flavors {uf = 2 + 1) of light quarks using domain- wall valence quarks on asqtad- 
improved staggered sea quarks [HI [T2], and with two flavors (n/ = 2) of light quarks using 
twisted-mass quarks [T3] and improved Wilson quarks [llHITj. These determinations are 
in agreement with the Roy equation values. The first calculation of the 7r+7r+ scattering 
phase shift was carried out by the CP-PACS collaboration, who exploited the finite-volume 
strategy to study s-wave scattering with n/ = 2 improved Wilson fermions [HI [15] at pion 
masses in the range 171^^ ~ 500 — 1100 MeV. The amplitudes obtained from the Lattice 
QCD calculations were extrapolated to the physical mass using a polynomial dependence 
upon the pion mass, instead of using the known pion-mass dependence of the amplitude 
based upon the symmetries of QCD encapsulated in xPT. In a recent paper, the Hadron 
Spectrum Collaboration (HSC) studied the s-wave 7i'^7t^ phase shift with pion masses in 
the range m^r — 390 — 520 MeV [T8]. Further, they have provided the first Lattice QCD 
calculation of the 7i~^7i~^ phase shift in the d-wave (/ = 2) [18]. 

In this work, which is a continuation of our high statistics Lattice QCD explorations [19]- 
123] , we determine the 7r+7r+ scattering amphtude below the inelastic threshold. Calculations 
are performed with four ensembles oi Uf = 2 + 1 anisotropic clover gauge-field configura- 
tions at a single pion mass of m.,^ ~ 390 MeV with a spatial lattice spacing of bg ~ 0.123 fm, 
an anisotropy of ,^ ~ 3.5, and with cubic spatial volumes of extent L ~ 2.0, 2.5, 3.0 and 
3.9 fm. Predictions are made for a number of threshold parameters which encode the lead- 
ing momentum-dependence of the scattering amplitude, and dictate the scattering length, 
effective range and shape parameters in the effective range expansion (ERE) of the inverse 
scattering amplitude. The Lattice QCD predictions are found to be in agreement with the 
Roy-equation determinations of the threshold parameters and phase shift, and with the 
available experimental data. Beyond the threshold region, the LEC's that contribute to 
the two-flavor chiral expansion of the scattering amplitude are determined, allowing for a 
prediction of the phase-shift at the physical pion mass to be performed at next-to-leading 



order (NLO). The predicted phase-shift is in agreement with the experimental data. 

The Maiani- Testa theorem demonstrates that S-matrix elements cannot be determined 
from stochastic lattice calculations of n-point Green's functions at infinite volume, except 
at kinematic thresholds |21] • Liischer showed that by computing the energy levels of two- 
particle states in the finite-volume lattice, the 2 — )■ 2 scattering amplitude can be recov- 
ered [25H3^ . These energy levels are found to deviate from those of two non-interacting 
particles by an amount that depends on the scattering amplitude (evaluated at that energy) 
and varies inversely with the lattice spatial volume in asymptotically large volumes. In 
this paper, Liischer's method is used to extract the phase shift from the lattice-determined 
energy-levels. 

This paper is organized as follows. In Sec. [II], we provide some details of the lattice 
calculations: we discuss the anisotropic clover lattices that are used and the determination 



of the anisotropy parameter. Sec. Ill gives a summary of the eigenvalue equation which is 
relevant to extracting phase shifts from lattice-measured energy levels, in the center-of-mass 
(CoM) system and in boosted (lattice = "laboratory" ) systems. The results of the Lattice 



QCD calculations are presented in Sec. |IV] and relevant fits that are used to determine the 
effective range parameters, up to and including the shape parameter, are discussed. Sec. |V] 
includes a summary of the relevant xPT formulas, the chiral fits to the lattice data, and the 
prediction for the n'^ir^ phase shift up to the inelastic threshold at the physical pion mass. 
Finally, a summary of our predictions and a discussion of the systematic uncertainties is 



given in Sec. VI 



II. DETAILS OF THE LATTICE QCD CALCULATIONS 



A. Anisotropic Clover Lattices 

Anisotropic gauge-field configurations have proven useful for the study of hadronic spec- 
troscopy [35H38] . and, as the calculations required for studying multi-hadron systems rely 
heavily on spectroscopy, we have put considerable effort into calculations using ensembles 
of gauge fields with clover-improved Wilson fermion actions with anisotropic lattice spacing 
that have been generated by the HSC In particular, the nj = 2 + 1 flavor anisotropic clover 
Wilson action [5S1 HD] with stout-link smearing [H] of the spatial gauge fields in the fermion 
action with a smearing weight of p = 0.14 and Up = 2 has been used. The gauge fields 
entering the fermion action are not smeared in the time direction, thus preserving the ultra- 
locality of the action in the time direction. Further, a tree-level tadpole-improved Symanzik 
gauge action without a 1 x 2 rectangle in the time direction is used. 

The present calculations are performed on four ensembles of gauge-field configurations 
with L^ X T of 16^ x 128, 20=^ x 128, 24^ x 128 and 32^ x 256 lattice sites, with a renormahzed 
anisotropy ^ = bg/bt where 6^ and bt are the spatial and temporal lattice spacings, respec- 
tively. The spatial lattice spacing of each ensemble is bg = 0.1227 ± 0.0008 fm |37] giving 
spatial lattice extents of L ~ 2.0,2.5,3.0 and 3.9 fm respectively. The same input light- 
quark mass parameters, btmi = —0.0840 and btnis = —0.0743, are used in the production of 
each ensemble, giving a pion mass of m,r ~ 390 MeV. The relevant quantities to assign to 
each ensemble that determine the impact of the finite lattice volume and temporal extent 
are m.,^L and ttIt^T, which are given in Table |T| In addition, we tabulate the pion masses 



on the four lattice volumes. As discussed in detail in Ref. [22], exponential finite- volume 
corrections to the pion masses are negligible for these volumes, a necessary condition for the 
application of Liischer's finite-volume method for obtaining phase shifts. Additionally, the 
predicted exponential finite- volume corrections to vrvr scattering near threshold are expected 
to be negligible [32] • Multiple light-quark propagators were calculated on each configuration 
in the four ensembles. The source locations were chosen randomly in an effort to minimize 
correlations among propagators. 



TABLE I: Results from the Lattice QCD calculations in the four lattice volumes, t.l.u denotes 
temporal lattice units. 



L^ xT 


16^ X 128 


20^ X 128 


24^ X 128 


32^ X 256 


L (fm) 


~2.0 


~2.5 


~3.0 


~3.9 


niT^L 


3.888(20)(01) 


4.8552(84) (35) 


5.799(16)(04) 


7.7347(74) (91) 


m^T 


8.89(16)(01) 


8.878(54)(22) 


8.836(85)(02) 


17.679(59)(73) 


m^ (t.l.u.) 


0.06943(36) (0) 


0.06936(12)(0) 


0.06903(19)(0) 


0.069060(66)(81) 



B. Determination of the Anisotropy Parameter, ^ 

In the continuum and in infinite- volume, the energy-momentum relation for the pion is that 



of special relativity, E 



mz 



IpI . In Lattice QCD calculations, this relation is more 



complicated due to the finite lattice spacing (including the violation of Lorentz invariance) 
and the finite- volume, resulting in E"^ being a non-trivial function of p, which has a polyno- 
mial expansion at small momentum. Retaining the leading terms in the energy-momentum 
relation, including the lattice anisotropy C,, the energy and mass in temporal lattice units, 
and the momentum in spatial lattice units (s.l.u) are related by 



-btE^i\n\)f = {hm^f + 



1 f2Tih, 



e 



L 



n 



(1) 



The Lattice QCD calculations of the energy of the single pion state at a given momentum 
p = ^n (where n is an integer triplet) allows for a determination of ^, and hence establishes 
the single-particle energy-momentum relation that is crucial for determining the scattering 
amplitude from the location of two-particle energy eigenvalues. We obtain ^ = 3.469(11) 
where the statistical and systematic uncertainties have been combined in quadrature. This 
is consistent with the value determined by Dudek et al. of ^ = 3.459(4) [18]. A fit to a 
higher order polynomial provides a result that is consistent with this value but with larger 
uncertainties in the contributing terms. It is important to use the lattice determined value of 
^, and to propagate its associated uncertainty, as small variations in this parameter are am- 
plified in the determination of the scattering amplitude from two-particle energy-eigenvalues 
when the interaction is weak (and the energy of the two-particle state is consequently near 
that of the non- interacting system) . 



III. THE FINITE VOLUME METHODOLOGY 



The formalism that was put in place by Liischer to extract two-particle scattering amplitudes 



below the inelastic threshold from the energy-eigenvalues of two-particle systems at rest in 
a finite cubic volume [271 I2H] was extended to systems with non-zero total momentum by 
Rummukainen and Gottlieb [30] . Subsequent derivations have verified and extended [ISl EU^ 
the work in that paper. The use of boosted systems allows for the amplitude to be 



determined at more values of momentum (in the CoM), between those defined by ^n. Here 
the results that are relevant to the present analysis of the boosted tt+tt^ systems, and to 
systems at rest, are restated. 

Using the notation of Ref. [31], the energy in the CoM frame is denoted by E*, which is 
related to the energy E and momentum Pcm in the "laboratory system" (the total lattice 
momentum) by E*"^ = E^ — |Pcmp. In what follows, it is useful to define Pcm = I Pcm I- The 
7-factor is straightforwardly defined by 7 = E/E*, and E* is also related to the magnitude 
of the momentum of each 7r+ in the CoM frame q* by E*^ = 4[ q*"^ -|- m^ ]. The real 
part of the inverse of the s-wave scattering amplitude below inelastic threshold, and hence 
the scattering phase-shift, can be extracted from the total energy of the two-particle system 
with total momentum Pcm = ^d in the finite- volume via the generalized Liischer eigenvalue 
relation 



q* cot S{q* 



'jLy/n 



^oo(i;?~ 



;;*2\ 



(2) 



where the dimensionless quantity q* is defined by q* = ^q* 
generalization of the functions defined by Liischer [271 128] , 



The function Zqq( 



■ ;7*2n 
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^LM(i;r') = E 
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g*2 
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where the Ylm are spherical harmonics and the sum is over vectors defined by 



1 f l^^ --u IjN 

- ( ny --d 1 +nx = 7 (n--d) 



(4) 



which in turn are related to the lattice momentum-vectors by k = ^n = ^ (ny + n^). The 
n are triplets of integers and the decomposition of n is along the direction defined by the 
boost-vector d. Liischer presented a method [271 l2H] which can be used [30] to accelerate 
the numerical evaluation of the sum in eq. (|3|, and a generalization of that method leads to 



ZtAM;q 



;;*2\ 



E 



-A(|r|2-g*2) 



1^ ^ /-/^ Zi 



|r|^ YLui^r) 
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dt (-j e*"? e ^ 



(5) 



where 



7W = 7W|| + w_L 



(6) 



The value of the sum is independent of the choice of A, and A 
works [TSl. 



1 has been used in previous 



The energy-level structure resulting from Z^Q(l;g*^) has been discussed previously, e.g. 
Ref. [30]. For the present calculations of boosted systems it is important to identify the 
closely spaced energy-levels. This is because the amplitudes extracted from such levels 
are subject to large systematic and statistical uncertainties due to the rapid variation of 
Zqq{1; q*"^) in their vicinity, and also due to the difficulty in separating the states contributing 
to the correlation functions. The energy-levels associated with two non-interacting particles 
are located at the poles of Zqq{1; g*^), and eq. Q gives 

d = (0,0,0) : q*' = 0, 1, 2, 3, 4, 5, .... 

H = m n n ~*2 ^ ±_ 47^ + 1 _9_ 87^ + 1 47^ + 9 167^ + 1 

'^ ' ' ^ ■ ^ 4^2' 4y ' 4^2' 4^2 ' 4^2 ' 4^2 ' •••■ 





472 ' 472' 


472 ' 

^v ^ 


472 


1 


272 + 1 3 


47^ + 1 


4 + 72 


2 


' 272 ' 2' 

^ V ' 


272 ' 


272 



d = (0,1,1) : r = ^, ^, ^^;^, :^, ^^, ^7^ •- , (7) 

and so forth, where the underbraces denote states that become degenerate as 7 — )■ 1. We 
stress that the relations summarized in this section are only valid below inelastic threshold. 

IV. 7r+7r+ SCATTERING ON THE LATTICE 
A. Lattice Phase Shift 



The scattering of pions in the 1 = 2 channel is perturbative at low momentum and at small 
light-quark masses, as guaranteed by xPT. In a finite volume, this translates into two-pion 
energies that deviate only slightly from the non- interacting energies; i.e., the sum of the pion 
masses (or boosted pion masses for moving systems). We have analyzed 7r+7r+ correlation 
functions with Pcm = 0, 1, a/2 and with various (non-interacting) momentum projections 
among the pions. It is straightforward to partially diagonalize this system of correlation 
functions into the energy-eigenstates at intermediate and long times. This is achieved by 
assuming that the two-pion energy levels are close to their non-interacting values, and then 
varying the linear combination of correlation functions in order to maximize the plateau 
region. (Coupled exponential fits to the various correlators with the same Pcm lead to 
consistent determinations.) As an example, in fig. [I] we show the two-pion effective mass 
plots (EMP's) on the 32^ x 256 ensemble with Pcm = 0, 1, \/2. Six energy levels can be 
clearly identified in the EMP's in fig. [I] for each of the values of Pcm- (Note that these levels 
clearly show the near degeneracies of the non-interacting system as established in eq. (^.) 
However, only the first few levels, when propagated through the eigenvalue equation, lead 
to statistically significant values for the phase shift. While the energies of other levels 
are established, the structure of the eigenvalue equation is such that the uncertainties, as 
small as they appear, are sufficiently large to produce uncertainties in the amplitude that 
are too large and preclude statistical significance^. The states that have been analyzed to 
produce amplitudes and phase-shifts are given in Table [Tl| and are shown in fig. [2j Note 

^ The EMPs of fig [I] indicate that the signal-to-noise ratio of the two-pion correlation functions decreases 
with increasing excitation number. 
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FIG. 1: The two-pion EMP's for the first six levels (here n indicates the level) with Pcm = (top), 
1 (middle) and v2 (bottom) in units of the temporal lattice spacing on the 32^ x 256 ensemble. 
Only one half of the temporal lattice points are shown. 



that momenta are quoted in units of rrivr in order to formulate the subsequent analysis in a 
manner that is independent of the scale setting. The values of k cot 5/m^ and 6 resulting 

TABLE II: Results from the Lattice QCD calculations of tt^tt^ scattering in the four lattice 
volumes. Pcm denotes the magnitude of the momentum of the center-of-mass in units of 2'k/L. In 
the column denoted by "level", g.s. denotes the ground state, 1st denotes the first excited state 
and 2nd denotes the second excited state. 



k^/ml 


L^^T 


Jem 


level 


k cot 5/m^ 


5 


0.00678(54)(81) 


323 ^ 256 





n = 


-4.49(35)(52) 


-1.06(12)(18) 


0.01772(14)(23) 


243 X 128 





n = 


-4.24(32)(49) 


-1.82(19)(30) 


0.0309(17)(27) 


203 X 128 





n = 


-4.25(21)(34) 


-2.37(18)(29) 


0.0715(32)(48) 


163 X 128 





n = 


-3.80(15)(22) 


-4.03(25)(35) 


0.1641(20)(23) 


323 X 256 


1 


n = 


-3.33(38)(48) 


-7.1(0.8)(1.0) 


0.378(5)(11) 


203 X 128 


1 


n = 


-4.1(0.4)(1.0) 


-8.6(0.8)(3.6) 


0.3838(42) (85) 


323 X 256 


^/2 


n = 


-1.65(12)(28) 


-20.6(1. 5)(3.1) 


0.7323(53)(88) 


323 ^ 256 





n = 1 


-2.78(29)(57) 


-17.2(1. 7)(2.9) 


0.9233(51)(73) 


323 X 256 


1 


n = 1 


-2.14(16)(26) 


-24.1(1.6)(2.6) 


L373(13)(22) 


243 X 128 





n = 1 


-2.10(19)(36) 


-29.2(2.3)(4.3) 


1.582(9)(16) 


323 X 256 





n = 2 


-1.19(09)(14) 


-46.5(2.3)(3.5) 


1.969(02)(04) 


203 X 128 





n = 1 


-2.33(32)(56) 


-31.6(3.5)(5.6) 



from the energy-eigenvalues are shown in fig. |3} 

Note that while the 323 ^ 256 P^m = \/2, ri = and 203 x 128 Pcm = 0, ra = 1 levels appear 
discrepant, we believe this is likely a statistical fiuctuation. Also, the phase shift we have 
extracted from the first excited state in the 243 ^ ]^28 ensemble disagrees with the equivalent 
extraction presented in Ref. p^8j. While we find a phase shift oi 5 = —29.2 ± 2.3 ± 4.3° at 

A;2 ~ 0.21 GeV^ Ref. [E] finds 5 13 ± 2° at k'^ ~ 0.2 GeVl Our result is consistent 

with the phase shifts at the nearby momenta calculated on the 323 ^ 256 ensemble. 



B. The Effective Range Expansion Parameters 

The ERE is an expansion of the real part of the inverse scattering amplitude in powers of 
the CoM energy. 



k cot (5 



1 



m. 



ra^^a 



m-jrV 



mt 



P {m.„r) 



k' 



mt 



+ 



(8) 



where m.,^a and mT^r are the scattering length and effective range in units of m'^, and P 
is the shape parameter. Here k = |k| is the magnitude of each pion's momentum in the 
CoM. Such an expansion is expected to be convergent for energies below the t-channel cut, 
which is set by vrvr exchange in the t-channel. The t-channel cut starts at fc^ = ?ti^, while 
the inelastic threshold is A;^ = Sm"^. 

As the calculations of k cot (5/?7i^ are approximately linear in k"^ in the region k'^/m'^ < 0.5, 
the scattering length and the effective range are fit (Fit A) using eq. ([s]) with P and the 
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FIG. 2: The two-pion energies in units of the temporal lattice spacing for the lattice ensembles 
considered in this work. The (vertical) thickness of each level indicates the uncertainty of the 
energy determination. Each state is labeled according to its center-of-mass momentum Pcm, and 
its excitation level n. The non-interacting levels are denoted by dashed (black) lines. Notice that 
the 32'^ X 256 Pcm = V^, n = and Pcm = 0, n = 1 levels are nearly degenerate. 



other higher order terms set to zero. The extracted values of ttIttO and ?Ti^r are given in 
Table III, and the resulting fit is shown in fig. |5| along with the 68% confidence interval 
error ellipses for the two-parameters. In the region k'^/m'^ < 1 the Lattice QCD calculations 
exhibit curvature consistent with quadratic (and higher) dependence on k"^. In Fit B the 
three leading ERE parameters are fit to the results of the Lattice QCD calculations. The fits 
are compared to the Lattice QCD calculations in fig. [6| which also shows the 68% confidence 
interval error ellipse for the two-parameter subspace of the three-parameter fit. It is clear 
from Table |III| that the fit parameters are consistent within the combined statistical and 
systematic uncertainties. In what follows, where we use xPT to predict the parameters at 
the physical point, the spread in value of the ERE parameters will serve as a useful gauge 
of the systematic uncertainty introduced in the fitting of the scattering amplitude. It is 
noteworthy that the data allows a significant determination of the shape parameter, P. 



TABLE III: ERE parameters extracted from the Lattice QCD calculations of kcotS/m-,^. 



Quantity 


Fit A: k'^/ml < 0.5 


Fit B: k'^/ml < 1 


rriT^a 


0.230(10)(16) 


0.226(10)(16) 


m^r 


12.9(1.5)(2.9) 


18.1(2.4)(4.7) 


rn^ar 


2.95(20)(42) 


4.06(30)(57) 


P 


— 


-0.00123(30)(55) 


xVdof 


0.83 


0.79 



10 





-1 

-2 
o -3 
-4- 
-5 
-6 



- 








- 


i K 






* 


.. 


1 \ 






■ ■ 


i 


► 


] 


I 






^ 


r 


1 


1 






• L/b,= 16 


- 










m L/b,= 20 


^ 


1 








A L/6,= 24 


- 


1 1 1 1 1 1 1 1 






♦ L/b,= 32 


1 


0.0 




0.5 


1.0 


1.5 


2.0 



FIG. 3: Results of the Lattice QCD calculations processed through the energy-eigenvalue relation 
to give values of kcoiS/rriT^. The {circles, squares, triangles, diamonds} ({black, red, blue, green}) 
correspond to the ensembles {l&^ x 128, 20^ x 128, 24^ x 128, 32^ x 256}. Statistical and systematic 
uncertainties are shown as the inner and outer error-bars, respectively. The vertical (blue) line 
at k'^ = ml 
inelastic threshold is at k'^ = 3m?, 



indicates the limit of the range of validity of the ERE set by the t-channel cut. The 



V. CHIRAL INTERPOLATIONS 



Motivation 



Although these Lattice QCD calculations have been performed only at one value of the pion 
mass, as we will see, the effective range and threshold scattering parameters satisfy low- 
energy theorems mandated by chiral symmetry, and therefore each scattering parameter 
can be used to fix the corresponding LEG that appears at NLO in xPT. Thus the scattering 
parameters at the physical point can be predicted at NLO in xPT. This is, in a sense, a 
chiral interpolation rather than an extrapolation since one is interpolating between the pion 
mass of the Lattice QCD calculation and the chiral limit. Unfortunately, the pion decay 
constant, f.„, has not yet been accurately computed on the anisotropic lattice ensembles that 
have been used in this work. However, xPT and the results of mixed-action Lattice QCD 
calculations |13] can be used to determine /^ (and its uncertainty) evaluated at the pion 
mass of the present L attice QCD calculations up to lattice spacing artifacts. Specifically, in 
what follows we use Vz^ = m^f^ff^ = 2.59(13) at m^ ~ 390 MeV. 
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FIG. 4: Results of the Lattice QCD calculations processed through the energy-eigenvalue relation 
to give values of the phase-shift 5. The phase-shift at low energies is shown as an inset. The {circles, 
squares, triangles, diamonds} ({black, red, blue, green}) correspond to the ensembles {16'^ x 128, 
20^ X 128, 24^ X 128, 32^ x 256}. Statistical and systematic uncertainties are shown as the inner 
and outer error-bars, respectively. The vertical (blue) line at k"^ = m^ indicates the limit of the 
range of validity of the ERE set by the t-channel cut. The inelastic threshold is at A;^ = 3m^. 

B. Threshold Parameters in xPT 

The relation between the tt^tt^ s-wave scattering amplitude t{s) (= t£zo(s)) and the phase 
shift 5 is given by [1] 



tis) 



4 



1/2 



2i 



{e 



2iS{s) 



1} 



(9) 



where s = 4(1 + k'^ /m^) and k = |k| is the magnitude of the three-momentum of each tt^ in 
the CoM frame. The NLO scattering amphtude can be expressed in terms of three LEC's, 
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FIG. 5: Left panel: a two-parameter fit to kcotd/m-,^ over the region k'^/m'^ < 0.5 (Fit A). The 
{circles, squares, triangles, diamonds} ({black, red, blue, green}) correspond to the {16^ x 128, 
20^ X 128, 24^ X 128, 32'^ x 256} ensembles. The shaded bands correspond to statistical (inner-yellow) 
and statistical and systematic added in quadrature (outer-pink). Right panel: 68% confidence 
interval error ellipses in the {m.,^a)~^ - ^rriT^r space. The inner-solid ellipse and the outer-dashed 
ellipse correspond to the statistical and to the statistical and systematic uncertainties added in 
quadrature, respectively. 
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Ci, C2, andC4 [21 E]: 

m = 






31 



+ 
+ 



P 



1 m. 



3847r3 
2 ^ 301 



47r /2 V 11527r3 



14 
457r3 



19 



Ci 



5127r' 



1 ^ 7 



-C2 + 2167rC4 



ra. 



+ 



+ 



47r3/4 V 32 
1 /I 



m. 



^„.p , 5^. ^ ,„„ ^ ^^^ 



m^ + m^/c + k 



A;2 



167r3/4 V4 - 

87r3/4 \^ 16 - ^ 9 " 18 



A;2 + m^ 




/,2 _j_ ^2 
T^ log 



1287r3/4 V 12 P ' ^ 



' fc^+mj 



fc2 



'fc2_^ ^ 1 



(10) 



The Ci can be expressed in terms of the li = ["^{fi = /tt), the famihar low-energy constants 
of two-flavor xPT [2J , 



Ci 



1 

'2^ 



(4/1 + 4/2 + /3-/4) 



1 



1287r3 ' 

C2 = 327i{12h + 4I2 + 7 h - 3U) + ^ ; 

ovr 



a = 



^ (212/1 + 40/2 + 123/3-69/4) + "^^ 



51847r2 



6220807r4 



(11) 



The behavior of the amplitude near threshold {k"^ — > 0) can be written as a power-series 
expansion in the CoM energy 

Ret{k) = -m^a + k'^b + k^c + 0{k^), (12) 

where the threshold parameters b and c are referred to as slope parameters. Matching the 



threshold expansion in eq. (12) to the ERE in eq. ([8]) gives |1]: 

1 



m.„r 



2m'ib 

—— + 2m^a 



P 



{niT.aY \{m.„ay' — A{mT,aY + 8(m7ra)^ — 4 (m,rO + 2(m7ra)^) b m\ — 8 (6^ + m^ra c) m^] 



(mTra — 2(m7ra)3 + 26 m 



2^3 



(14) 
These equations can be inverted to obtain b and c from the lattice-determined ERE pa- 



rameters. Expanding the NLO amplitude in eq. (10) in powers of /c, one finds NLO xPT 
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expressions for the ERE and threshold parameters: 

z 2^ 3;z2 247r ^ 17, 

m^a = ^ + 2; Ci + ———log 2 , m^r = + C2 + -^logz , 

2 o ^ 11^ , r. 23^' 3^ 613^3 

"^-"^ = 3 + zC3 + ^log. , P = -^^^^^ + z C4 + ^^^^^^log^ , 

/ 10 q 5 

m^c = -.= (-C. - ^^G + 216.C. + ^log. ) . (15) 

where z = m^/ f"^ and C3 = 247rCi + g^C2. It is important to note that the shape parameter 
P and the threshold parameter c do not receive contributions from LO xPT; i.e. they vanish 
in current algebra. 



C. Chiral Interpolation of Threshold Parameters 



Using the ERE parameter set from Fit B given in table III , with statistical and systematic un- 
certainties combined in quadrature, the four functions Ci((m7r'2)'"**5 -z'""), C2((Tn^r)'° **, 2: '°**), 
C3((m^ar)'°", z'°"), C4(P''*", 2;'"") can be determined. The ERE parameters in table III give 



fjNLO ^ -0.00237(52) , C^^^ = 5.2(5.2) ; 

Cg^^o = -0.02(0.10) , Cf^^ = 9.0(4.0) x 10"^ , (16) 

where the superscripts denote that these constants are evaluated at NLO in xPT, and from 
which follow, using eq. (15), the predictions at the physical point^ of 



m^a = 0.0417(07)(02)(16) , m^r = 72.0(5.3)(5.3)(2.7) , mlar = 2.96(11)(17)(11) , 
P = -2.022(58)(12)(76) x 10^^ , 
b = -0.832(50)(0)(31) x 10~^ m^^ , c = 0.013(33)(01)(0) m^^ , (17) 

where the first systematic uncertainty has been estimated by comparing the interpolated 
results of Fits A and B, and by "pruning" the highest energy datum from the Lattice QCD 
results and refitting. The second systematic uncertainty contains estimates, using naive 
dimensional analysis (NDA), of the effects from higher orders in the chiral expansion, next- 
to-next-to-leading order (NNLO) and higher, as well as the contributions from lattice spacing 
artifacts that are expected to contribute at 0{b1) |11]. The chiral interpolation of mlar is 
shown in fig. [TJ (Note that the band in fig. [T] represents Fit B, and the outer uncertainty 
on the interpolated result represents the effect of the two systematic uncertainties and the 
statistical uncertainty added in quadrature, as described above.) 



With eq. (11), the fit values of the C, in eq. (16) can be used to constrain various 



combinations of the /j, and the renormalization group can be used to express these constraints 



^ Note that the precise NPLQCD result for the scattering length, m,ra = —0.04330(42), computed in 
Ref. [T^ with domain- wall valence quarks on staggered sea quarks, is more precise than the result of 
eq. (|17|). 
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FIG. 7: The dashed (green) Hne denotes the physical line, and the horizontal solid (purple) line 
denotes the LO xPT prediction, which is rn^ar = 3 in the chiral limit. The band denotes the 68% 
confidence interval interpolation of the results of the lattice calculation (the (red) rectangle) from 
Fit B. The Lattice QCD +xPT prediction at the physical point is the (red) star on the physical 
line, and the Roy equation prediction [8j is the (black) circle on the physical line. 



in terms of the scale-independent dimensionless barred quantities, the /j [2] ^- We find that 



h - 4/4 
2I1 - 3/3 



-29(27) 
28(29) , 



, /i — 6/4 
I1 + 4I2 = 



= -32(25) 
15.8(6.7) , 



(18) 



where statistical and systematic uncertainties have been combined in quadrature. With 
increased precision in the determination of the ERE parameters, such determinations of the 
LECs could become competitive with other methods. 

These results may seem surprisingly accurate for a Lattice QCD calculation performed 
at a single pion mass. As mentioned previously, it is the chiral symmetry constraints on 
the scattering parameters in the approach to the chiral limit that is responsible for the 
precision. The scattering length obtained here is consistent within uncertainties with the 
previous Lattice QCD determinations pTV[T3] . Further, the scattering length and threshold 



■^ The L are related to the L via 



I, 



32tt 



h + log 



where 71 = i, 72 = |, 73 = 



— —^ and 74 = 2. 
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parameters are found to agree with determinations from the Roy equation (with chiral 
symmetry input) [8], 

m^a = 0.0444(10) , b = -0.803(12) x lO^^m^^ ; 
mlar = 2.666(0.083) , (19) 



at the la-leveL Fig. u\ provides a comparison of the lattice calculation (and interpolation) 
and the Roy equation value of mlar. 

D. Chiral Interpolation of the Phase Shift 

The TT^vr"'" scattering phase-shifts calculated with Lattice QCD, which extend above the 
range of validity of the ERE but remain below the inelastic threshold, can be used to 
predict the phase shift at the physical value of the pion mass. While the chiral expansion 
may break down for scattering at sufficiently high energies, we ignore this issue and fit the 
NLO xPT amplitude (one-loop level) to the results of the Lattice QCD calculations at all 
of the calculated energies, the maximum invariant mass being ^/s ~ 1340 MeV. 

The results of the Lattice QCD calculations given in Table |Tl] are fit to the formula 

kcotS _ I, , k"^ f 1 iNLo(^) ^ , • ^ fr,n,^ 



m^ \j ml \ ti,o{k) {t'Lo{k)Y J rrij, 

where ^lo and tNLO are the LO and NLO contributions to t{k) in the chiral expansion, given 



in eq. (10). The result of the fit is shown in fig. [8] in the left panel the fit (of Ci, C2, 
and C4) to kcot6/m,T, is shown, and in the right panel, the fit values of Ci, C2, and C4 
(fully correlated) are used to predict the phase shift at the pion mass of the Lattice QCD 
calculations, ?ti^ ~ 390 MeV. The 68% confidence intervals for Ci, C2, and C4 from this fit 
are 

^NLO ^ (-0.0040,-0.0013) , C^^^ = (2.67,24.1) , 

fjNLO ^ (_i 7^_^3Q) X 10"^ , (21) 

with a x^/dof = 2.1 (for the fit with the statistical and systematic uncertainties combined 
in quadrature). The interpolated ERE parameters are: 

m^a = 0.0412(08)(16) , m^r = 80.0(9.58)(3.0) , P = -1.85(31)(07) x 10~^ , (22) 

which are consistent within uncertanties, but less precise, than the threshold determinations 



of eq. (17). Here the second uncertainty is an NDA estimate of the effects of higher orders in 
the chiral expansion and lattice spacing artifacts. For a better determination of the threshold 
parameters from the global fit, one requires more accurate Lattice QCD calculations and the 
7r+7r+ amplitude beyond NLO in the chiral expansion. In fig. [9] the fit values of Ci, C2, and 
C4 are used to predict the phase shift at the physical value of the pion mass, m^r ~ 140 MeV, 



which is compared to the experimental data of Refs. [45H48J . Fig. 10 compares the phase 
shift prediction to the Lattice QCD phase-shift determination by CP-PACS [U], and the 
Roy equation determinations of the phase shift from Refs. [El [9] . One should keep in mind 
that the interpolated phase shift is valid above the inelastic threshold, as the An intermediate 
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state appears beyond NLO in the xPT calculation (at two-loop level). The combined Lattice 
QCD and xPT prediction of the vr+vr"'" phase shift at the physical pion mass is found to be in 
good agreement with the experimentally-determined phase shift. While for |k|^ 400 MeV 
the uncertainty in the prediction exceeds the uncertainties in the experimental data, below 
this momentum the Lattice QCD+^PT prediction is more precise. 
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FIG. 8: Three parameter fit (Ci, C2, and C4) of the NLO xPT expression for k cot 5 /ra-,^ to the 
results of the Lattice QCD calculations. The shaded bands correspond to statistical (inner-yellow) 
uncertainties and statistical and systematic uncertainties added in quadrature (outer-pink). The 
solid (black) curve in the right panel is the LO xPT prediction (current algebra) at the pion mass 
used in the Lattice QCD calculations, m-,^ ~ 390 MeV. 

It is interesting to observe that while the LO phase-shift well reproduces the results of 
the Lattice QCD calculations at m,r ~ 390 MeV, as shown in fig. |8| the NLO contributions 
are important at the physical pion mass, as seen in fig. 



10 This is to be contrasted with 



the chiral behavior of the scattering length which is dominated by the LO amplitude, with 
NLO making a small but noticeable contribution. In an attempt to isolate the origin of this 
apparent difference, it is useful to consider scattering at NLO in the chiral limit where 



cot 5 — )■ 



4vr5 
A;2 



+ 



40 , f2k 
9vr \f\ 



9 224 

SSvr^Ci - — C2 + MhQ-K^Ci - -— . (23) 

62 AoTT 



The phase-shift can be defined this way even in the chiral limit because at LO and NLO 
the only intermediate states contributing to the scattering amplitude involve two pions. 
Inelastic channels, such as four-pion intermediate states which would invalidate the relation 
in eq. (lo]), first contribute to the scattering amplitude at NNLO. This is what allows for 
the phase-shift to be predicted above the inelastic threshold, and to remain perturbatively 
close to the actual value for momenta below the chiral symmetry breaking scale. At LO 
in the expansion, the phase-shift reaches 6 = 7r/4 when k'^ = 47r/^ ~ 0.22 GeV^ 
Jtt = 132 MeV), consistent with the phase-shift shown in fig 
to take the limit m^ ^ k for this value oi k {k r^ 470 MeV). 
k, the NLO terms are approximately equal to the LO terms, providing an estimate of the 
convergence region of the chiral expansion for the scattering process. 

It is also worth noting that while it is formally invalid to use the Liischer relation in eq. ^ 
for the scattering of pions above inelastic threshold, xPT indicates that the error introduced 
into phase-shift determinations is small, occurring at NNLO in the chiral expansion. This 
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FIG. 9: The shaded band is the Lattice QCD prediction of the phase shift at the physical value of 
the pion mass, m-,^ ~ 140 MeV using NLO xPT with the statistical and systematic uncertainties 
combined in quadrature. The data is experimental (black and grey) taken from Refs. |45H48| . The 
red vertical line denotes the inelastic (47r) threshold. 



is not expected to be true for other scattering processes (those not involving the pseudo- 
Goldstone bosons). Therefore, while strictly speaking the results presented in Ref. [18] 
above inelastic threshold arise from an invalid application of eq.([2]), the expected deviation 
from the true result is expected to be small (at momenta for which the chiral expansion is 
converging), suppressed by two orders in the chiral expansion. Clearly, precision calculations 
of the phase-shift above the inelastic threshold cannot rely upon a methodology that does 
not include the effects of inelastic processes. As all of the calculations in our work are 
below the inelastic threshold, the present analyses and predictions do not suffer from this 
inconsistency. 



VI. SUMMARY AND CONCLUSION 

The increases in high-performance computing capabilities and the advent of powerful new 
algorithms have thrust Lattice QCD into a new era where the interactions among hadrons 
can be computed with controlled systematic uncertainties. While calculation of the basic 
properties of nuclei and hypernuclei is now a goal within reach, it is important to consider 
the simplest hadronic scattering processes as a basic test of the lattice methodology for 
extracting scattering information (including bound states) from the eigenstates of the QCD 
Hamiltonian in a finite volume. In this work, we have calculated the tt^tt^ scattering 
amplitude using Lattice QCD over a range of momenta below the inelastic threshold. Our 
predictions for the threshold scattering parameters, and hence the leading three terms in 
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FIG. 10: The shaded band is the Lattice QCD prediction of the phase shift at the physical value 
of the pion mass, m.,^ ~ 140 MeV using NLO xPT with the statistical and systematic uncertainties 
combined in quadrature. The outer (magenta) band is the CP-PACS physical prediction |15j , and 
the two inner (blue and purple) bands are the Roy equation predictions [U |9]. (Note that the 
band from Ref. [8j lies above the band from Ref. |^.) The solid (purple) curve is the LO xPT 
prediction (current algebra) at the physical pion mass. The red vertical line denotes the inelastic 
(47r) threshold. 

the ERE expansion, are consistent with determinations using the Roy equations [H |9] and 
the predictions of xPT. In particular, our determination of m^ar = 2.96(11) (17) (11) from 
an interpolation of a fit to the low momentum values of k cot 6/171^^ is consistent with the LO 
prediction of xPT of m^ar = 3 (l + Oim^/ h?)). Further, the resulting predictions for the 
phase shift at the physical pion mass -using NLO xPT- are in agreement with experimental 
data, and are even more precise in the low-momentum region. 

The Lattice QCD calculations presented here were performed at one lattice spacing sim- 
ply due to the lack of computational resources, therefore, an extrapolation of the ERE 
parameters to the continuum limit (as was performed in the work of CP-PACS p^) could 
not be performed. The discretization of the quark fields that has been employed gives rise 
to lattice spacing artifacts at Oipl), and we expect such contributions to be small for these 
calculations. 
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